Relative velocity of sliding of microtubules by the action of Kinesin-5 
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Kinesin-5, also known as Eg5 in vertebrates is a processive motor with 4 heads, which moves 
on filamentous tracks called microtubules. The basic function of Kinesin-5 is to slide apart two 
anti-parallel microtubules by simultaneously walking on both the microtubules We develop an 
analytical expression for the steady-state relative velocity of this sliding in terms of the rates of 
attachments and detachments of motor heads with the ATPase sites on the microtubules. We first 
analyse the motion of one pair of motor heads on one microtubule and then couple it to the motion 
of the other pair of motor heads of the same motor on the second microtubule to get the relative 
velocity of sliding. 



I. INTRODUCTION 

The protein constituents of the cytoskeleton of eu- 
karyotic cells can be divided broadly into the follow- 
ing three categories: (i) filamentous proteins, (ii) motor 
proteins, and accessory proteins. The three classes 
of filamentous proteins which form the main scaffold- 
ing of the cytoskeleton, are: (a)ocim, (h) microtubule, 
and (c) intermediate filaments. The three superfamilies 
of motor proteins are: {i)myosin superfamily, iii)kinesin 
superfamily, and {lii) dynein superfamily. Microtubules 
serve as tracks for kinesins and dyneins.pj. 

The main cellular mission of the molecular motor 
Kinesin-5 is to cross-link anti-parallel microtubules and 
to slide them apart, thus playing a critical role during 
bipolar spindle formation j^j . Kinesin-5 is absolutely es- 
sential for the dynamic assembly and function of the mi- 
totic and meiotic spindles f3|. 

In this work, I have closely followed the paper by Pe- 
skin and Oster 



II. MODEL 

The molecular motor Kinesin-5 is a four headed motor 
with two heads on each side. FIG-[T] shows a schematic 
representation of the model we have developed for this 
motor. In our model, we assume the two heads of the 
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FIG. 1. Schematic diagram of Kinesin-5 on microtubules 

motor on one side to be joined at a hinge. Similarly, 
the other two heads on the other side are also joined 



by another hinge. These two hinges are then linked by a 
flexible connector. We also assume that the two heads on 
each side are free to rotate about their respective hinge. 
The two heads can also pass each other freely. However 
there is one restriction on their relative displacement. 
The maximum separation between two heads on one side 
can be equal to the separation between two discrete bind- 
ing sites on the microtubule. This separation, which is 
a constant, is denoted by d. We assume that the micro- 
tubules are aligned anti-parallel because Kinesin-5 has a 
much stronger preferences for cross-linking anti-parallel 
overlapped microtubules than parallel ones (3]. The -I- 
end directed stepping activity of Kinesin-5 j^] slides mi- 
crotubules away from each other. FIG E] shows schemat- 
ically the model of one pair of heads on a microtubule. 
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FIG. 2. Schematic diagram showing only one side of the motor 
with its two heads at discrete adjacent binding sites on a 
microtubule. 

The factors included in the model which are responsi- 
ble for the motor to move forward are: (i) The back head 
has a higher probability to detach than the front head, 
(ii) The attached head prefers a geometry in which it 
leans in the forward direction i.e towards the -|- end of 
the microtubule. The first of these factors ensures that 
the back head detaches and the second one ensures that 
it reattaches as a front head f5|. When both the heads 
are attached to the microtubules then the back head is in 
a more favourable orientation than the front one. There- 
fore, ATP hydrolysis and subsequent detachment is more 
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likely at the back head. Once the back head detaches, the 
front head relieves itself of the strain by leaning forward 
and swinging the detached back head forward. Hence this 
is the process by which the motor moves forward. The 
coordinates of the two heads are taken to be a: = xi and 
X — X2- The hinge is located halfway betwween them. 
The coordinate of the hinge, xu is therefore 



Xh 



Xl + X2 



and the restriction dictates |xi — a;2| < d 
Let 

a — rate constant for attachment of the free head 
to the track 

rate constant for detachment of the back head 
from the track 

rate constant for detachment of the front head 
from the track 
= probability that the free head binds infront of 

the bound head 
- p = probability that the free head binds behind 
the bound head. 

Each head is also assigned a number which signifies its 
state. If it is attached to the track it is assigned 1 and if 
it is detached from the track it is assigned 0. 

The state 00 corresponds to both the heads of the cor- 
responding side completely detached from the track and 
that implies the termination of the walk. So we only con- 
sider the states 11, 01 and 10. FIG13] shows the possible 
transitions and their rates. 
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FIG. 3. Possible transitions in the states of one pair of heads 
of the motor and the transition rates 

In order to monitor the motion of one pair of heads 
on a microtubule and the velocity of this motion, it is 
simpler and adequate to monitor the motion of the mid- 
point of the heads Q . The average velocity with which 
the mid-point moves on the microtubule can be taken as 
the average velocity with which the microtubule slides 
past the motor heads. We define, 

Xm = coordinate of the mid-point, which is identified 
as follows [5|: 

Xm = position of the hinge in the 11 state 

Xm = position of he bound head in the 10 = 01 

state 



So when both the heads are attached then the mid- 
point is halfway between the two ATPase sites. As soon 
as one head detaches, the mid-point shifts to the location 
of the bound head i.e one of the two binding sites on 
either side of the previous location of the mid-point. So 
from FIGlHwe can see that the midpoint jumps ±| and 
this is the step size of the motion of the mid-point. A 
state diagram showing the possible steps of the mid-point 
and their rates is shown in FIG|4l It is assumed that the 
binding sites are labelled by half-integer indices and the 
mid-points between them are, consequently, labelled by 
integer indices. So when both the heads are attached, the 
mid-point is at an integer site and when only one head is 
attached, the mid-point is at a half-integer site. 
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FIG. 4. The steps the mid-point can take and their respective 
rates 

So using the idea of Peskin and Oster explanined 
above, the motion of a pair of heads on one side can 
be simplified to the motion of the mid-point of the head 
alone. Similarly, the same thing can be done for the 
other pair of heads as well. So instead of monitoring the 
motion of four heads we can monitor the motion of the 
two mid-points like a two-headed motor. Now these two 
mid-points are coupled together because they are parts of 
the same 4-headed Kinsesin. This constraint introduces 
a restriction on the difference between the coordinates of 
the mid-points. If Xmi and Xm2 are the coordinates of 
the mid-points then restriction is \xmi — Xm2\ ^ f • So 
this ensures that if the mid-point is at a site i on a mi- 
crotubule then cross-linking can take place only through 
sites i and z ± i on the other microtubule. Similarly, if 
the mid-point on a microtubule is at a site i + 5, then 
the possible cross-linking sites on the other microtubule 
are i + ■^ji and i + The orientations allowed are shown 
in FIG 



III. CALCULATIONS AND RESULTS 

Let Pj {t) = probabilty of finding the mid-point at j at 
time t. From FIGH] we can write the following master 
equations: 

dPj{t) 



dt 



apPj_i + a{l — p)P, 



if3f+MPj (1) 
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Mean position = ^a;P(a;) 
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FIG. 5. The possible orientations of cross-links. The squares 
represent discrete sites where mid-points can be present. The 
circles represent the mid-points. 



dt 



/?/P,+i+/3fcP, -aP^-+i (2) 



Now we define 

Mo = E,^-oo and No = E,^-oo Pj+^- 



dMo _ ^ dPj 
dt ~ ~dt 



^ dt 

— — OO 

oo 



dt 



= aNo - (/?/ +;3fc)Mo 



(3) 



Similarly, 



dNo 
dt 



dt 



{(3b + /3f)Mo-aNo 



(4) 



In the steady state ^ = ^ ^ Q. So from equations 
© and gl) we get 

aiVo = + pb)Mo (5) 
Also from normalisation we have, 

Mo + iVo = 1 (6) 
Solving equations ^ and (H]), we get 



Mo 
No 



(7) 
(8) 



MeanPositic 



oo -j^ 

E bPj + ij + i;)Pj+i. 



MeanSpeed 



1 oo 

= ^[Eb-^. + (^- + 2^^.>^]] 



If we define < w > as the average speed of the mid- 
point then 



<v> 



using the above relation and equations dD) and © we 
can write 



< V > 



oo 



j=-oo 



Using Mo = Ej°l-oo and iVo = Ej°l-co ^j+i' get 



Mo 



. Plugging in the expressions for Mq and A'o from equa- 
tions and dH), we get (Sj] 



<v > 
d 

where, 



a p 



a + j 



2{a + pf + l3b) 



(9) 



l^Pb + Pf 
S = (3b~l3f 

Now we couple the motion of the two mid-points which 
move on their MT tracks with < v > calculated above. 
The motion of a mid-point counts only if it is attached to 
a microtubule and cross-linked to the other microtubule. 
So the relative velocity due to the motion of a mid-point 
on a microtubule is equal to < w > multiplied to the 
probability of a mid-point being found at any particular 
site on that microtubule, which is then multiplied by the 
probability of a mid-point being found at all the corre- 
sponding possible cross-linking sites, and then summed 
over all the sites of the first microtubule. 

This prescription can be mathematically expressed as, 
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using the definitions of equations (1121) and ([13] 



<v^ >=<.>, + + 

i— — oo 
i— — oo 

+ (10) 

Here < u >i and < v >2 are the velocities of the mid- 
points on the first and second microtubules respectively. 
Similarly, the superscripts on the probabilities also refer 
to the the microtubules. But, since there is no functional 
difference between the two microtubules, these indices 
can be dropped. Then we have 



<V>^2<v>[Y, [Pl+Pl^+^P, 

i— — oo 

We define three new quantities, 

oo 

i—~oo 
oo 

^= E 

i — — OD 

OG 

Q- E p^ 

i— — oo 

Plugging these expressions in equation ([TT|) . 

<V>^2<v> {AM + N + Q) 

Now, 

dN 



(11) 

(12) 
(13) 
(14) 

(15) 



2 E P^H 



dt ^ dt 

i—~oo 



using equation ([2]) 
dN 



dt 



dN 
~dt 



^2[{pf + Pb)M -aN] 



(16) 



Similarly, 



rfQ _ ^ p dPr 

dt 2^ dt 

2 — — oo 



using equation ([T]) 



dQ 



^=2 ^ [P.(apP,_i +«(1 - p)P,+ !-(/?/ +/3.)P0] 

2 — — oo 

using the definitions of equations ([T2|) and (fT4)) 



^=2[aM + ail-p)M-{pf+pb)Q] (17) 



In the steady state ^ = ^ = 0, so we have, 



2[{l3f + l3b)M-aN]^0 (18) 
2[aM + a{l-p)M~{f3f + f3b)Q]^0 (19) 



Also since the variables M,N and Q capture all the cross- 
linking probabilities (FIGE)) so we can write. 



2M + N + Q = 1 



Solving equations (UHl), (0 and ^ we get. 



(20) 



M 



2a(/3/ + /3b) + i/3f+/3b)^ + a^2~ p) 



(21) 



Using equations ([20| and ((2T|). we get the final expression 
for < F > as. 



<y> = 2<i>>(l + 2M) 



I.e., <V > ^2 



ad 

a + 7 



7 P 



2a{l3f+M 



2a{l3f + pb) -f (/?/ + /36)2 + a2 (2 - p) 



(22) 
(23) 



IV. PHYSICAL INTERPRETATIONS 

From equation (PT|) it is obvious that Af > 0. There- 
fore, <V>>2<v> which is also consistent with our 



intuitive expectation. Since in our model cross-linking is 
allowed through sites next to the adjacent sites on mi- 
crotubules, hence the two mid-points move with respect 
to each other. This factor actually increases the relative 
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velocity above 2 < v >. The microtubules slide with re- 
spect to the mid-points with veocity < v > and the mid- 
points move with respect to each other which increases 
the overall relative velocity. 

Next we consider a few special cases to establish that 
our results are consitent with physical intuition. 

If the probability of the detachment of the back head 
and the front head is the same and the probability of the 
free head binding infront and behind the attached head 
is also the same, then 



P =h 



(24) 
(25) 



In that case, using equations ([24l) and ([25]) in equations 
dH) and ([23]), we get < v >= and < V >= respec- 
tively, which is expected on physical grounds. 

If possible orientations for cross-linking are restricted 
by imposing the condition that cross-linking is allowed 
only through adjacent sites on two microtubules, then 
the relative positions of the mid-points with respect to 
each other do not change. The relative velocity of slid- 
ing will simply be the sum of the velocities with which 
each mid-point hops on the microtubule. The velocites 
are added because of the anti-parallel orientations of the 
microtubules. Mathematically, in this case 



M = £ P,P,+ i = 



(26) 



and hence, using equation (|26l) in equation ([2^ we get 

<V">=2<w>, (27) 
which agrees with the above prediction. 

V. COMPARISON WITH EARLIER WORKS 

Kruse et al. Q calculated the relative velocity of sliding 
between two filaments due to the action of a two headed 
motor. In this work the length of the filaments, which is 
denoted by L is assumed to be finite. It is assumed that 
cross-linking can take place only through adjacent sites 
of the two microtubules. A variable ^ is defined, which is 
the differnce between the coordinates of the minus ends 
of the microtubules. Thus overlaps can take palce if ^ = 
0, 2L — 2. The average velocity of sliding as a function 
of ^ in units of monomers per unit time is 

. / Y.\ti J^S2+i-^ for^ = 0, i - 2, 

v[^) - 'P^^ I Et"/-L+2 J^S2+,-. fore = L - 1, 2L - 3. 

(28) 

where pd is the probability of cross-linking. 



When averaged over all possible values of ^ the final 
expression for average relative velocity is Q 



2L-3 



<V>= 



2L 



(29) 



Now, equation can be used for our system by taking 
the limit of L ^ oo and replacing v{^) hy 2 < v >. In 
that case, 



2L-3 



2L-3 



lim — 

L— foo 2Ij 



<v> 



2 < V > ^^^^ 
lim -TT 1 

2L-3 



2L 



= 2 < V > lim 

L— >oo Z_L — 1 

<V > ^2 <v > 

which is consistent with equation (j23|) because in the spe- 
cial case cross-linking was allowed only through adjacent 
sites on two microtubules. 

Computer simulations as well as analytical studies re- 
lated to continuum models of motor induced sliding of mi- 
crotubules and actin filaments have been reported Q . 
These studies say that two headed bipolar motors slide 
anti-parallel filaments with a velocity that is twice that 
of the free motor velocity Q . In our case, the motion of 
Kinesin-5 is reduced to a motion of a two headed motor 
using the model of mid-points. In our case cross-linking 
is allowed through adjacent sites on microtubules or sites 
next to adjacent sites on microtubules. In the continuum 
limit the distance between two sites is negligible, so it can 
be said that all the cross-linking that takes place is effec- 
tively through adjacent sites only, as in the continuum 
limit the adjacent site and the next to adjacent site can 
be deemed to be the same. Hence all the cross-linking 
can be deemed to be through adjacent sites. In such a 
case as we calculated in equation ((^7)) our realtive veloc- 
ity of sliding < > is equal to 2 < u > , where < u > is 
the velocity with which a mid-point slides a filament. If 
< u > is taken analogous to the free motor velocity of a 
two headed motor, then our result in the continuum limit 
agrees with earlier studies done with continuum models. 
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